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In the Yang-Mills theory, the apparent thickness of the confining string is known to grow loga- 
rithmically when its length increases. The same logarithmic broadening also happens to strings in 
other quantum field theories and domain walls in statistical physics models. Even in quantum field 
theories, the observables used to measure and characterise this phenomenon are largely borrowed 
from statistical physics. In this paper, we describe it using the string form factor, which is a mean- 
ingful quantum observable, and show how the form factor can be obtained from field correlation 
functions calculated in lattice Monte Carlo simulations. We apply this method to 2+1-dimensional 
scalar theory in the strong coupling limit, where it is equivalent to the 3D Ising model, and through 
duality also to 2+1-dimensional Z2 gauge theory. We measure the string form factor by simulating 
the Ising model, and demonstrate that it displays the same logarithmic broadening as statistical 
physics observables. 



PACS numbers: 05.50.+q, 11.27. +d 



I. INTRODUCTION 

String-like excitations play an important role in many 
quantum field theories, for example the confining string 
in Yang-Mills theory [l[ and cosmic strings in some Grand 
Unified Theories String excitations are also closely 
related to the physics of interfaces in statistical physics @, 
1. 

In the semiclassical approximation, the string is de- 
scribed by a solution of the field equations which is time- 
independent and translation invariant along the string, 
and describes the properties of the string on all length 
scales. Quantum mechanically, the picture is very differ- 
ent because the string carries massless Goldstone modes 
whose quantum fluctuations dominate the dynamics on 
small scales. For example, the string width appears to 
depend logarithmically on its length because of these fluc- 
tuations. 

Because the behaviour of the Goldstone modes is inde- 
pendent of microscopic details, this logarithmic broaden- 
ing is a general property of any strings irrespective of the 
theory, and it has been studied extensively in Yang-Mills 
theory an d sp in models, both analytically @-0|, and nu- 
merically [10H16| . In the context of the confining string in 
QCD or Yang-Mills theory, it is generally seen as a prob- 
lem, because it makes it impossible to measure "intrinsic" 
properties of the string which one would like to know for 
phcnomenological descriptions of confinement 0, [13] • 

The previous numerical studies of strings in quantum 
field theory have been generally based on approaches that 
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are natural in statistical physics models, such as measur- 
ing expectation values of the plaquette near the string. 
The observables available in a quantum field theoretical 
system are generally quite different. If we really had ac- 
cess to a real, physical string, we would probably probe 
it experimentally by scattering particles off it. This mo- 
tivates us to focus on the string form factor, which is 
related to the corresponding scattering amplitude and is 
a well-defined quantum observable. Building on previ- 
ous work on kinks and monopoles [HI, [l9| , we show how 
the form factor of a string-like topological soliton can be 
calculated in Monte Carlo simulations. 

Because quantum mechanical strings in Euclidean 
spacctime are equivalent to domain walls in statistical 
mechanics, we apply the method in practice to calcu- 
late the domain wall form factor in the three-dimensional 
Ising model near the critical point. This theory is in the 
same universality class as the 2+1-dimensional real scalar 
field theory and has therefore the same critical behaviour. 
It is also exactly dual to the confining three-dimensional 
Z2 gauge theory, so our conclusions should also be valid 
for confining strings, at least qualitatively. 



II. STRING SOLUTION IN SCALAR THEORY 

Let us start by considering the 2+1-dimensional real 
scalar field theory with the Lagrangian 
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where /1 e {0,1,2}. In the broken phase, which semi- 
classically corresponds to m 2 < 0, the theory has two 
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vacua 

* = ±v = ±^. (2) 
Classically, there is also a topologically stable solution 
4>{t, x, y) = v tanh ^ m (^_^) ^ ^ ^ 

which we have chosen to be perpendicular to the y direc- 
tion. In 2+1 dimensions, this solution can be thought of 
either as a string (because it is a one-dimensional object 
on any time slice) or a domain wall (because it divides 
spacetime into two pieces). Since we are principally in- 
terested in using these objects to model strings in 3+1- 
dimensional theories, we will refer to it as a string when 
we discuss it in the context of the 2+1-dimensional the- 
ory. 

There is a general result Q, applicable to any string- 
like objects, that fluctuations broaden the string so that, 
whatever its initial shape, its width w is given by 



where a is tension of the string. 

However, the arguments used to demonstrate this are 
based on the dependence on the expectation value of the 
energy density on the transverse position relative to the 
string, or other similar quantities which are sensitive to 
the fluctuations of the string position. There has, there- 
fore, been discussion on whether there is some other ob- 
servable which would measure the "intrinsic width" of 
the string HQ. 

III. STRING FORM FACTOR 
A. Definition 

A stable quantum state corresponding to Eq. ([3]) also 
exists in the quantum theory. We can define it as the 
lowest topologically non-trivial energy eigenstate, and we 
denote it by |0). To be precise, there is a degenerate set 
of such states corresponding to different orientations of 
the string, and we again choose the string that is oriented 
along the x direction as shown in Figure [TJ 

The string ground state |0) has zero momentum and, 
in the infinite- volume limit, its energy is Eq = aL, where 
the constant a is the string tension. We can obtain mov- 
ing string states by boosting the ground state |0) in the 
y direction. We denote these states by \p), where p is 
the momentum of the string in y direction. The energy 
of such a state is E p = \J p 1 + Eq. We normalise these 
states as 

(p'\p) =2n6(p'- P ). (5) 

This normalisation is not Lorentz invariant under rela- 
tivistic boosts in the y direction, but because the string 
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FIG. 1. System setup. The coordinate x runs along the do- 
main wall (shaded gray) on a timeslice, the twisted boundary 
conditions are in the y direction and t parameterises the di- 
rection along which correlators are measured. 

has a linearly divergent energy, we can safely restrict our- 
selves to only non-relativistic motion. 

For any local operator O, such as the field operator </>, 
we can now define the corresponding form factor as 

f(p2,Pi) = (p 2 |O(0)bi), (6) 

where by O(0) we mean the operator in coordinate space 
at position (x,y) = (0,0). 



B. Semiclassical Limit 

In the semiclassical limit J2(j, the form factor is given 
by the Fourier transform of the classical profile O c \(y) of 
the quantity O. Taking matrix elements of the Heisen- 
berg equation of motion for the classical field we find 

/(pa,Pi) = (P2|O(0)|pi) 

= j dy e^-^y<D ci (y) 

= O c \{p 2 - Pi)- 

(7) 

This expression is valid only in the non-relativistic limit, 
when \pi\, \p%\ *C <jL. In this limit the form factor is a 
function of the momentum difference k = p2 — pi only, 
as a direct consequence of the Galilean invariance, so wc 
will denote it by f(k). Although our studies of kinks and 
monopolcs used fully rclativistic semiclassical expressions 
for the form factor, this is a reasonable approximation in 
the current paper since any results obtained with a string 
moving relativistically are most probably due to finite- 
size effects (the amount of energy required to accelerate 
a string to relativistic velocities on macroscopic scales is 
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too great; the finite box size does not correctly capture 
this). 

In particular, choosing O — cj>, wc find from Eq. © 
the semiclassical result for the string form factor 



/(*) 



(8) 



m sinh(fc7r/\/2m.) 
At low momenta k -C m, this behaves asymptotically as 



/(*)=«*, 



(9) 



which corresponds to an infinitcsimally thin string 

My) = uSign(y). (10) 

Because any string would look like this when viewed from 
long distances, we expect the low-fc asymptotic behaviour 
(j9|) to be valid for any string, irrespective of microscopic 
details and even in quantum theory. Therefore it serves 
as a useful benchmark for our calculations. 

C. Correlator and form factor 

It was shown previously in Ref. [HI how the form fac- 
tors of point-like solitons such as kinks and monopolcs 
can be calculated from the field correlation function. The 
same technique works also for strings. 

Wc work in momentum space, taking the Fourier trans- 
form in space but not in the time direction. 

After Wick rotation to Euclidean time, the two-point 
correlator of the operator O has the spectral expansion. 



(0(Q,k x ,k y )0(t,q x ,q y )) 

(0\6(k x ,k y )\a}{a\d(q x ,q y )\0) 



-t{E a -E ) 



(0|0) 



(11) 



where U(t) = exp(-Ht) is the Euclidean time evolution 
operator, and Z — Tr U(T) the partition function for 
the worldsheet, not the full Zt w of the field theory with 
twisted boundary conditions. In other words, only when 
T — > oo is the situation of Eq. (fTTj) realised; otherwise 
the string's worldline is restricted. 

With twisted boundary conditions, the states \a), \a') 
must have an odd number of strings (in practice, exactly 
one due to heavy volume suppression factors e~ 2aLT 
in the partition function; this also suppresses opposite- 
spin bubble formation). Because of momentum conser- 
vation they must also have opposite overall momentum 
(k x ,k y ) = —{q x ,q y )- If the momentum is in the y di- 
rection, i.e., k x = 0, then the lowest such state is the 
moving string state \k v ) with momentum k y , as defined 
in Section fill Al It has energy 
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where Eq is the ground state energy of the string. 

In contrast, if the momentum has a non-zero x compo- 
nent, k x 7^ 0, then the interpolating states must be ex- 
cited states with some excitation carrying the momentum 
in the x direction. The lightest such states are massless 
Goldstonc modes which exist because the position of the 
string breaks translation invariance in the y direction. 
Because we are only interested in the unexcited states 
\k), we do not consider this possibility, and instead we 
simply restrict ourselves to k x = from now on. To sim- 
plify our notation, we therefore suppress the unneeded 
argument k x , and write 



O(ky) = d(0,ky), 



(14) 



and 



where Eq is the energy of the single string ground state. 

A small added complication is that the spacetime is 
necessarily finite in actual Monte Carlo simulations. We 
have periodic boundary conditions in the ^-direction (as 
well as the x-direction along the string), and twisted 
boundary conditions in the y-direction. This periodic- 
ity leads us to write the correlator as 

(O(0,k x ,k y )O(t,q x ,q y )) 

= -i-Tr U(T - t)6(q x ,q y )U(t)6(k x ,k y ) 



^^2(a'\d(q x ,q y )\a}(a\d(k x ,ky)\a 



l\ e -E a ,(T-t)-E a t 



(12) 



(O(0,ky)O(t,q y )) = (O(0,0,ky)O(t,0,q v )). (15) 



The next states in the spectrum contain pairs of mass- 
less Goldstone modes, with opposite quantised momenta 
k x = 2rm/L. These propagate along the string per- 
pendicular to our chosen direction for time evolution of 
the system. The energy of the lowest such state, with 
k x = 2tt/L, is Ek y +4ir/L. Therefore, if we choose t > L, 
these states are strongly suppressed relative to the mov- 
ing string state in the spectral expansion (|12| . Conse- 
quently, the problem beco mes p ractically identical to the 
l+l-dimensional kink case [18| . and we can approximate 
Eq. (|T2|) by an integral over moving string states \k y ) 
only, 
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dk' dk" 

(O(0,k y )O(t,q y )) = - I ^^(k'\d(q v )\k'')(k''\d(k v )\k 



-E k ,(T-t)-E k 



2ir 2ir 



L 2 

— 2ir5{q y 



V V 



\2 e -E k , y (T-t)~E kly _ ky t 



(16) 



where we have used the result 
(k y \6(q v )\%) = / dxdye^y{k' y \d{x,y)\k^ 

dx dye iq y y e^ k 'y-K)y{k' y \d(0)\k^) 

= 2nL8(q y + k> y ~ki>)f(ki y ,kii). (17) 

The choice t > L obviously requires a relatively long 
lattice in the time direction. The dynamics of the 
longest-wavelength modes then becomes essentially one- 
dimensional. 

Similarly, we can write the partition function as 



dk 



Z= ^-{k'\U{T)\k')=L /-f 



L 



2- 

dk 



dk' 



-E k , T 
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This partition function is the individual contribution to 
the partition function from string's classical worldsheet. 

To calculate the integral (fir?]) , we use the saddle point 
approximation. The saddle point ko is found by minimis- 
ing the action 



S(k' y )=E K (T-t) + E k 



t 



(19) 



for given t. By approximating the integral by a Gaussian 
around the saddle point, we obtain 



(O(0,k y )O(t,q y )) = ^2Tr5(q y + k y )x 



dk' , , 
2ir y 



ky,k y )\ 



2 -S(to)-iS"(to)(fc;-fco) 2 



(20) 



This is not a Gaussian approximation of the correlation 
function, nor is it a scmiclassical stationary phase calcu- 
lation; it does, however, allow us to rearrange Eq. (|16|) 
and solve for \f(k' y , k' — k y )\ 2 in exchange for imposing 
minor restrictions on what we can measure. In the limit 
of large t and T — t that we already require, the Gaus- 
sian approaches a delta function and we can calculate the 
integral 



(O(0,k y )O(t,q y )) 
■■ 2iT5{k y +q y )L 



T 
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point fco by minimising Eq. (|19p , and we obtain 



f(k , fc - ky) = ±iy/{O(0,ky)O(t, 

'E S"(k Q 



-ky)) 



T 



1/4 



(S(k )-E a T)/2 



(22) 



for O odd. Note that the saddle point fco still depends 
on t. 

The result (f2"2")) should be compared with the corre- 
sponding result for kinks given in Eq. (19) of Ref. [lH 
(although we have corrected a typographical error here). 
In that case, the interesting length scales were compa- 
rable to the inverse kink mass. Therefore the motion of 
the kink was relativistic for the corresponding momenta, 
and it was natural to express the result in terms of ra- 
pidities. In the current case, the interesting momenta 
ky ~ \fa are much lower than the string mass uL, and 
its motion is therefore highly non-relativistic. We can 
therefore simplify Eq. (|22[) by taking the non-relativistic 
limit. In that case, fc = (t/T)k y and S"(k ) = T/E , 
and the form factor also becomes a function of the the 
momentum difference k y only, 

f(k y ) = /(fc , fco - k y ) 

' kl t(T-t)\ 



= ±U (O(0, k y )0{t, -ky)}- exp 



4£o 



T 



(23) 



Correspondingly Eq. (|2T|) simplifies to 
(O(0,k y )O(t,q y )) 



2 7 rL5(fc !/ + g y )|/(fc)| 2 cxp 



l- 2 

v 



t(T - t) 



2En 



T 



(24) 



We can also use this result to determine the ground state 
energy Eq from the correlator measurements. 



D. Linear fluctuations 



|/(fc , fc -fcy)|V 



We can use Eq. ([2T]) to determine the form factor from 
the field correlator. For given k y and t, we find the saddle 



While the primary purpose of this paper is to mea- 
sure the form factor fully non-perturbativcly using lat- 
tice Monte Carlo simulations, it is useful to look at the 
S(k )+E T leading-order quantum effects analytically. The lowest 
excitations in the system are massless Goldstone modes 
(2$i) the string worldsheet, which are present because the 
string breaks translation invariance in the y direction 
spontaneously. To calculate their effect, we assume a 
string that moves in the normal direction y without 
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changing its profile. If the profile of the string is </> c (y), 
and the position of the string is y(t,x), then the field 
configuration is 

<t>(t,x,y)=4>c(y-y(t,x)). (25) 
Taking the Fourier transform over x and y, we find 

Mt,k x ,k y ) = 4> c {k y ) [ dxe^ x+k ^^\ (26) 



where (j> c (k y ) is the Fourier transform of the classical pro- 
file, 



Mk v )= I dye ik ^4> c {y). 



(27) 

We can now write the correlation function as 

(0(t, k X ,kyW, k' x , k' y )) = ^ c {k y )^ c (k' y ) 

x f dxdx' e i( - k °> x+k '* x,) ^(k y y(t,x)+k' v y(t',x'))^ _ ^ 

Writing Ax = x' — x and At = t' — t, we have 

(4>{t,k x ,k y )(j>{t' ,k' x ,k'y)) = 2ir S(k x + k' x )4> c (k y )^ c (k y ) 
x J dAxe~ lKAx ^ e i(k y y(o,o)+k' y y(At,Ax))y (29) 

Translation invariance in the y direction gives a delta 
function also for the k y component, so we can assume 
k' y = —ky and we have 

(Cf>(t, k X ,kyW, k' x , ~ky)) = 2%8{k x + k' x )\4> c {k y )\ 2 

x J dAxe- zk * Ax ^(k y (y(0,0)+y(AtAx)))y (30) 

Assuming that y(t, x) is Gaussian, this is equal to 

(0(t, k x ,k y )(/)(t', k' x , -k y )) = 2n6(k x + K)\4> c (k y )\ 2 

X / d & x e -^Ax e -±kl((y(0,0)-y(At,Ax)) 2 ) , ( 31 ) 



However, if we consider an elongated lattice with spa- 
tial size L and time separation At L, then the cor- 
relator is very different. Let us start with the two-point 
correlator of the field y, which is given by the Fourier 
transform of the propagator 

, , . . „ 1 f dk f duj e l ( kx +^ , . 

(vMy(t, X )) = - j - J -j^, (32) 

and in finite spatial volume, the integral over k becomes 
a sum over n, so we have 

I JfL r dbJ pi(krtX+Ult) 

{ y M y( t ,x)) = - £ ^^~r (33) 



with k n = 2im/L. This is obviously IR divergent, but 
we can use it to write 

((y(0,0)-y(t,x)) 2 ) =2(j/(0,0) 2 )- 2(y(0,0)y(t,x)) (34) 

which yields 



1 



t + 2j2 



1-e 



-k„t 



cos(k n x) 



((y(0,0)-y(t,x)Y) = — 

(35) 

For t > L, the exponentials are all very small and we can 
approximate 



((y(0,0)-y(t,x)) 2 )^ 
1 



OO ^ 



n=l 



7r ' n 

n=l , 



(36) 



(37) 



This diverges in the UV, but a minimum length scale / 
provides a cutoff n < Ljl. Then for L ^ ^ one has 



n=l 

and we find (for k x = 0), 

(^(t,0,kyW,0,-ky)) 

= L 2 \4> c (k y )\ 2 e-; 



ln 7+7 , 



(38) 



-(ln(L/£)+ 7 ) -(fc3/2<rL)i 



(39) 



Inserting this into Eq. (f2"5| . we find the form factor 

f(k v ) = $ c (k y )e-&WW+r). (40) 

Comparing with Eq. ([7]), we can see that the effect of the 
fluctuations is to suppress the form factor at high mo- 
menta k y > yj2-Ka ln(L/£). For large L, this washes out 
any structure the classical string solution has at short 
distances, and using the asymptotic low-momentum be- 
haviour of the classical solution </> c (fc y ) ~ 2iv/k y , we find 



/(*) = f e" 



r (ln(L/*)+ 7 ) 



(41) 



If interpret, in line with Eq. ([7|), this as the Fourier 
transform of the quantum-corrected domain wall profile 
0cff (y), we find that in coordinate space the domain wall 
has broadened and has width 



t<>~ 



1 

2na 



(ln(L/£)+ 7 ), 



(42) 



in perfect agreement with Eq. ([4J. 



IV. ISING MODEL 



A. The model 



To demonstrate the use of Eq. (|2"Tj) we use it to calcu- 
late the domain wall form factor non-perturbatively near 
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the critical point. This is an interesting calculation be- 
cause the theory becomes strongly coupled and therefore 
perturbation theory is not valid. In practice, we do not 
actually simulate the scalar field theory but rather the 3D 
Ising model, which is known to be in the same universal- 
ity class and which will therefore give identical results 
near the critical point. From a computational point of 
view, the Ising model is much more convenient because 
highly efficient numerical algorithms are available to sim- 
ulate it. 

The worldsheet of the string corresponds to a domain 
wall in the 3D Ising model. In this section we will there- 
fore use the term domain wall, but it should be under- 
stood to refer to the same physical object as the term 
string elsewhere in the paper. The 3D Ising model is also 
exactly dual to the confining 2+1-dimensional Z2 gauge 
field theory. This duality maps the domain wall to the 
worldsheet of the confining string of the gauge theory, 
and therefore our calculation also describes the proper- 
ties of the confining string in this somewhat simple gauge 
theory. 

The logarithmic broadening ^ of the domain wall is 
well known also in the Ising model. In the rough phase 
(between /3r and /3c); the domain wall width w has a 
logarithmic divergence with domain wall length L Q of 
the form given in Eq. ((4]). 

The Ising model has Hamiltonian 

H = -J2 fcx,x's( x M x (43) 

<x,x') 

where x and x' refer to position vectors in the 3D Eu- 
clidean space, and (x, x') denotes that the sum runs over 
all nearest-neighbour links in the lattice. The spins s(x) 
take the values ± 1 . We can denote the size of each direc- 
tion L x , L y and T respectively, but for our own calcula- 
tions we shall always take L x = L y = L. It has partition 
function 

Z = E exp(-0H), (44) 

{s(x)=±l} 

and we impose periodic boundary conditions in all three 
directions, with fc XiX ' = 1 for translational invariancc. 
However, to create a domain wall, a twist is introduced 
so that fc Xi x' = — 1 in the y-dircction for one value of y. 

Various techniques exist to measure the free energy of 
such an object; almost all depend on measuring the ratio 

F = -ln^, (45) 

where Zq is the partition function on a fully periodic 
lattice and Z tw is the partition function with antipcriodic 
boundary conditions in the y-direction. 

The free energy of the domain wall can be obtained 
as a function of L x and T . For a cubic system with 
T = L x = L, a perturbative calculation yields [2l[ 

Fcubic = 0L 2 - X - In o + G - — + O ((aL 2 )- 2 ) , (46) 



where G « 0.29, and the domain wall tension a cor- 
responds directly to the string tension in the 2+1- 
dimensional scalar theory. 

B. Real-space width measurements 

Previous measurements of the width of walls and 
strings have worked in real space. Here we summarise 
one of the more successful approaches and note that later 
developments are generalisations of the same idea to dif- 
ferent systems or improvements in numerical technique. 

In Ref. [22j], the domain wall width was studied. The 
results were based on linear functions of local spin op- 
erators; in addition the i-direction was not privileged in 
this study. For these reasons, this calculation cannot be 
considered a direct analogy of a particle scattering off the 
domain wall. 

Note that the geometry used was rather different to 
ours. In addition, the discussion here follows our own 
conventions established in Section [J First the total mag- 
netisation for slices parallel to the domain wall was mea- 
sured, 

For each configuration, m(y) was calculated, and the 
minimum value yo of y obtained. Then the 'normalised 
magnetisation gradient' is 

p(y) = JjTf\™>(y-yo + l) -rn(y-y )\ (48) 

where 

M = jJ2 m (y) ( 49 ) 

y 

is the average magnetisation. For a very sharp and 
straight domain wall, p(y) will only be nonzero for one 
value of y. In reality, fluctuations will smear the result, 
even for individual configurations. There may also be 
bubbles of spin in the bulk, or (for small lattice sizes) 
other odd numbers of domain wall. The height h can be 
defined as 

* = £>(»), (50) 
y 

and then the variance in this expression gives us the mean 
squared width w 2 , 

™ 2 = P^y 1 (E pMvJ ) : ( 51 ) 

the angle brackets denote an ensemble average. Note that 
if the shifting procedure docs not introduce a systematic 
bias in (h) (say, the value of yo is chosen randomly when 
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there arc multiple degenerate minima of m(y)), then we 
have (h) = which we shall assume shortly when study- 
ing how this estimate of w 2 behaves. 

Anyway, using this data for various lattice sizes (for 
example, in Ref. [l3l | the system under consideration had 
L x = L t and L y up to 27, so the width either side of 
the domain wall was fixed) , a fit to the ansatz taking the 
same form as Eq. ([4]) can be made. 

We assume (h) = 0, as discussed. Now 



E 



y 

2LM 



^2(s{t, x, y) - s(t, x, y - 1)) 



(52) 



(53) 



p=0.226 
Nambu-Goto 



■■'T 



ISO 
L 



In other words, this measurement technique is 
weighted towards long-distance fluctuations of the do- 
main wall. It does not and cannot directly probe the 
physics of any intrinsic width, nor does it allow for the 
possibility of different behaviour being visible at short 
distances as such effects will be washed out. One way of 
evading this problem is the rescaling approach adopted 
in Ref. [H. 

The approach discussed above is equally applicable to 
the confining string in SU(iV) gauge theory, as was first 
carried out in Ref. [l4|. Of course, the endpoints of the 
string must be fixed as well so there is not even transla- 
tional invariancc in such a calculation. 



V. RESULTS 

We use the standard Ising single-cluster algorithm [23| , 
as well as Metropolis updates to obtain the results pre- 
sented here. To show the volume-dependence of the 
ground state energy we use (3 = 0.226; however the 
domain wall is so light here that finite-size effects for 
the form factor for higher k are severe. Therefore the 
form factor results arc presented for values of j3 that are 
slightly further into the rough phase. 



A. Stationary string energy 

In order to directly recover the string form factor from 
correlator measurements, we need an independent esti- 
mate of the string tension a. This is obtained from mea- 
surements of the the lowest momentum (k = ir/L) corre- 
lator in the presence of the string. 

The expression in Eq. (|24| is used, along with a contri- 
bution from bulk scalar particles; these have the standard 
cosh-type long distance behaviour and their contribution 
to the correlator can be easily accounted for. To improve 
fitting performance, the scalar mass is obtained from fits 
to correlators with periodic boundary conditions, leaving 
only the relative amplitudes of the two contributions and 



FIG. 2. Energy per unit length of the string measured using 
the correlator with f3 — 0.226. Units are a 2 ; the error bars 
are jackknife estimates. 

the surface tension undetermined. This approach is in 
line with our previous studies along the same lines. 

Note that at higher momenta, the bulk scalar contribu- 
tion ~ g-w+Pf decays much more rapidly with t than 

the string contribution ~ e 2E o , so discarding the first 
few points at either end of the lattice allows us to ignore 
this issue when measuring the form factor; Figure [3] also 
shows this effect. 

As expected, we see is evidence for the Luscher term in 
our string tension measurements. For a closed Nambu- 
Goto string, the ground state energy is 



ML) = *L^l-— 2 . (54) 

Unfortunately, we cannot rule out the massless world- 
sheet modes affecting the quality of the data for the 
smallest of our lattice sizes. The results are shown for 
(3 = 0.226 in Figure along with a fit to Eq. (|54| . 

Nonetheless, the quality of this fit gives us confidence 
that it is principally the ground state of the string that 
gets excited when we carry out measurements at large 
time separation; this is in line with the assumptions of 
Section UrCl 

B. String form factor 

The previous section presented evidence to support our 
approximation that the correlator method we employ in 
this paper picks up only the ground state of the string. 
We proceed now to measure the form factor and interpret 
this as a way of measuring the interactions of the string 
on various length scales. 

The string moves non-relativistically k <C <jL in almost 
all the cases studied, so while the expressions derived in 
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FIG. 3. Typical plots of the string form factor for /3 = 0.23, 
L = 64, as a function of distance for several k. The error 
bars here are bootstrap estimators for the random error in 
the correlator measurement; they are therefore correlated at 
different values of t. Only the results for 16 < t < 48 are 
independent of separation and it is results in that region that 
are used in the other plots presented later in this paper. 



FIG. 4. String form factor for /3 — 0.23 and various L, 
with fits given by the ansatz in Eq. (|55[) . assuming the short- 
distance behaviour takes the form <j)(k) — 2/k, implying that 
the string has no discernible intrinsic width. In this plot and 
those that follow, the systematic measurement error in a is 
not shown on the a;-axis. 



Section IHI CI allow for this, there is no significant de- 
pendence on the form factor with time separation; there 
are only short-distance effects that could, in principle, be 
removed. 

Let us first demonstrate the applicability of Eq. (|2lT) . 
Figure [3] shows the form of the processed form factor 
results (we also tested the technique with T = 4L and 
quantitatively identical results were obtained, albeit at 
considerably greater computational cost). It is assumed 
that the form factor is pure imaginary, we have no way 
of determining its phase. We cannot subtract the bulk 
scalar contribution; this is the principal source of devia- 
tions from the expected behaviour at short and interme- 
diate distances. 

Taking results at long distance from correlators similar 
to those shown in Figure [3l we arrive at Figure [4] Here, 
appropriately scaled form factor measurements are plot- 
ted as a function of the dimensionless combination k/aL. 
For comparison, the step-function result from Eq. © is 
shown. A thin string of this form would be independent 
of volume L, and it is clear that at long distances (small 
k), the behaviour of the form factors is independent of 
volume. 

While Figure 5] is scaled in such a way that the re- 
sults can be compared with the form factors of pointlikc 
solitons (as the x-axis plots total momentum over rest 
energy) and to demonstrate that the string moves non- 
relativistically in most of the results presented here, to 
see whether the string form factor were independent of 
volume a different scaling combination is required. With- 
out logarithmic broadening, we would expect that the 
data points in Figure [5] would collapse onto a single line 





L=32 






L=64 






L=128 






Z.=256 
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FIG. 5. As for for Figure [4] but scaled instead by ^fa. One 
would expect the combination k/^/a to parameterise the in- 
trinsic width of a string. The data do not collapse onto a 
single curve. 



for all k. This is clearly not the case, and we instead 
have dependence on volume. Irrespective of the scaling 
employed we see that for larger L, the form factor falls 
off faster with increasing fc, suggesting that the string 
indeed becomes broader at larger volumes. 

Finally, Figure [6] shows that because w 2 (as defined 
in Section |TT| actually depends logarithmically on the 
correlation length £, there is no way to collapse the data 
onto a single curve. 
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FIG. 6. As for Figure [5] but for fixed L = 48 and various /3. 
As can be seen, the data do not collapse onto a single line 
even for constant L due to the weak dependence of w 2 on the 
bulk correlation length. This supports our ansatz, Eq. (|55|l . 

C. Broadening 

The analytic approximation in Eq. (|41| suggests that 
the string form factor should decay as 

|/(fc)|«|0(fc)|e-^ V (55) 

where w 2 is the 'width' of the string due to fluctuations 
as derived in Section [ill D I and 4>(k) is the intrinsic form 
factor of the unbroadened string. Unfortunately, when 
the string is moving non-relativistically, it is difficult to 
distinguish between different choices of 4>(k) - all of which 
must fall off as 2iv/k at small k. For the time being we 
shall assume that the simplest model - that of a step 
function at short distances, given in Eq. (fT0|) - is valid 
and take 4>(k) oc 1/k. 

If the string had intrinsic width and took the form 
of a smooth kink then (j>(k) would be given by Eq. ([8]). 
However, for k/(aL) < 1 at fixed L the difference in 
the two curves is insignificant and therefore we favour 
the simpler fit arising from step-function model, even if 
it does not capture all the physics on the shorter length 
scales. We shall discuss the prospects for measuring an 
intrinsic width in the conclusions. 

It is not possible to establish conclusively based on 
Figure |4] that logarithmic broadening is taking place; in- 
deed the high-/c behaviour of Eq. © is rather similar. 
However, in Figures [5] and [6] we can see that even when 
correctly scaled the data do not collapse onto a single 
line for different j3 or L. 

Finally, the results of fitting Eq. (f55"j) to the data shown 
in Figures [4] or [5] are shown in Figure [7] (although addi- 
tional values of L are shown in the latter). There is a 
clear linear relationship between w 2 and logL, consis- 
tent with Equations (|U) and (f3"9")) . We must therefore 




L 



FIG. 7. The values of w 2 obtained when fitting Eq. (T55|) 
to the data shown in Figure 0] this shows the width of a 
string as a function of L for /3 = 0.23. A linear relationship is 
apparent and a fit of these points to Eq. (|42p is shown, yielding 
a = 0.01 91 ± 0.0003 - a result that is fairly close to but not 
in excellent agreement with the literature, perhaps due to 
unaccounted systematic effects. Nonetheless, the behaviour 
is exactly that of Eq. ([4]) . 

conclude that the behaviour observed is entirely due to 
the fluctuations in the string. 



VI. CONCLUSIONS 

We have demonstrated how to study the properties of 
topological strings and domain walls in quantum field 
theory using two-point field correlation functions. In the 
actual numerical calculation we used the 3D Ising model 
which is in the same universality class as the 2+ ID real 
scalar field theory, and therefore has the same behaviour 
near the critical point. The model is also exactly dual to 
the confining Z2 gauge theory. 

The measured ground state energy Eq shows finite- 
size behaviour which is compatible with the Liischer term 
generated by quantum fluctuations of Goldstone modes. 
In the string form factor, these same Goldstone modes 
suppress short-wavelength modes leading to behaviour 
that can be interpreted as logarithmic broadening of the 
effective string width. This is consistent with previous 
studies of other measures of string widths or domain wall 
widths in statistical physics models. However, we stress 
that in contrast with the quantities studied before, the 
string form factor is a well-defined quantum observable. 

Although, in principle, it would be possible to use the 
form factor to probe the intrinsic structure of the string, 
this would require excellent data and may not be feasible. 
Any deviations from zero intrinsic width in the results 
presented here were within the errors of our fits.. 

For zero-dimensional topological solitons there are no 
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massless fluctuations and the spectrum of Goldstone 
modes is relatively straightforward. For example, in the 
case of a A</> 4 kink, the only Goldstone mode is that which 
we exploit in measuring the kink's mass. In the present 
model, this situation is complicated considerably. How- 
ever we have successfully demonstrated that our tech- 
nique is again of use. 

We would expect our approach to be applicable, with 
some modification, to confining strings in SU(iV) gauge 
theory. 
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